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ABSTRACT 

We present a new method to find voids in cosmological simulations based on the tidal 
and the velocity shear tensors definitions of the cosmic web. We use the fractional 
anisotropy (FA) computed from the eigenvalues of each web scheme as a void tracer. 
We identify voids using a watershed transform based on the local minima of the FA 
field without making any assumption on the shape or structure of the voids. We 
test the method on the Bolshoi simulation and report on the abundance and radial 
averaged profiles for the density, velocity and fractional anisotropy. We find that voids 
in the velocity shear web are smaller than voids in the tidal web, with a particular 
overabundance of very small voids in the inner region of filaments/sheets. We classify 
voids as subcompensated/overcompansated depending on the absence/presence of an 
overdense matter ridge in their density profile, finding that close to 65% and 35% 
of the total population are classified into each category, respectively. Finally, we find 
evidence for the existence of universal profiles from the radially averaged profiles for 
density, velocity and fractional anisotropy. This requires that the radial coordinate is 
normalized to the effective radius of each void. Put together, all these results show 
that the FA is a reliable tracer for voids, which can be used in complementarity to 
other existing methods and tracers. 

Key words: Cosmology: theory - large-scale structure of Universe - Methods: data 
analysis - numerical - N-body simulations 


1 INTRODUCTION 

Cosmic voids are regarded as one of the most striking fea¬ 
tures of the Universe on its larger scales ever since they 


von Benda-Beckmann & Muller 2008; 

Foster & Nelson 2009; 

Ceccarelli et al.|2013 

Paz et al. 2013 

Sutter et al. 2014 

1 . 


were found in the first galaxy surveys (Chincarini & Rood 
1975JJ&egory&Thompson[J1978JJEinastoe^^L 1980ajb[ 


On the theoretical side, the basic framework that ex¬ 
plains the origin of voids was established in the seminal work 


of Zel’dovich (1970) and refined in the following decades. 


Kirshner et al. |1981| Zeldovich et al. 1982| Kirshner et al.| The g rs t detailed theoretical models describing formation, 


1987). However, due to the large volume extension of void 


Redshift Survey (2dF) (Colless et al. 2001 20031 and the 
Sloan Digital Sky Survey (SPSS) {York et al.||2000 Abaza- 
|jian et al. 1 12003] >. These observational breakthroughs gen¬ 
erated a great interest in the last decade to study voids 


regions 5 — 10 Mpc/i ), statistically meaningful cata- 

Icke 1984 Bertschinger 1985 Blumenthal et al. 19921 were 

logues of voids (Pan et al. 2012 Sutter et al. 2012 Nadathur 

complemented and extended by numerical studies (Martel & 

& Hotchkiss 2014) have only become available through mod- 

Wasserman 1990, Regos & Geller||1991; |van de Weygaert & 


(Hoyle & Vogeley|2004| Croton & et al. 2004 

Padilla et al. 

2005 

|Rojas et al. 2005; Ceccarelli et al. 2006 

Patiri et al. 

2006 

Tikhonov 2006; Patiri et al. 2006; Ti 

honov 2007; 


Currently, the most popular approach to study voids relies 
on N-body simulations. For an extensive compilation of pre¬ 
vious numerical works we refer the reader to|Colberg et al. | 
{2008]). 

The relevance of voids to cosmological studies can be 


summarized in three aspects (Platen et al. 2007). Firstly, 
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voids are a key ingredient of the Cosmic Web. They domi¬ 
nate the volume distribution at large scales and compensate 
overdense structures in the total matter budget. Secondly, 
voids provide a valuable resource to estimate cosmological 
parameters as their structure and dynamics are sensitive to 
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them. Finally, they are a largely pristine environment to test 
galaxy evolution. 

Although visual recognition of voids in galaxy surveys 
and simulations is possible in most cases, we need a clear al¬ 
gorithmic identification procedure to make statistical stud¬ 
ies. Nevertheless, the community has not reached yet an un¬ 
ambiguous definition of cosmic voids. There are many differ¬ 
ent void finding techniques in the literature (for a detailed 
comparison of different schemes, see the publication on the 
results of the Void Finder Comparison Project|Colberg et al.| 


I2008)). In spite of the diversity of existing schemes, they can 


be roughly classified into two types: point-based and field- 
based. There are geometric schemes based on point distribu¬ 
tions (either real or redshift space) f galaxies in surveys or 


Muller et al.:2000. |Gottlo 

ber et al. 2003 Hoyle & Vogeley 

2004| 

Brunino et al. 2007 

Foster & Nelson 2009: Micheletti 

et al. 

2014 

Sutter et al. 

20151. While other schemes are 


based on the smooth and continuous matter density field ei¬ 
ther from simulations or from reconstruction procedures on 
surveys (Plionis fe Basilakos|2002' Colberg et al.|2005| Shan 


darin et al.|2006T Platen_et^aL|200T |Neyrinck|2008 JA unoz- 

Cuartas et al.||2011| Neyrinck et al.|[2013[ Ricciardelli et al. 


2013). Our work follows the tradition of the second kind of 


schemes. 

Here we introduce a new algorithm to define voids over 
the continuous matter density or velocity distribution de¬ 
fined on a fixed and homogeneous spatial grid. The algo¬ 
rithm uses the results from two tensorial schemes used to 
classify the cosmic web. The first (the T-Web) is based on 
the Hessian of the gravitational potential or tidal tensor 
( Hahn et al.||2007 Forero-Romero et al.||2009| ). The second 


(the V-web) is based on the velocity shear tensor (Hoffman 
et al. ] 2012 ). Our procedure allows for a description of a void 


internal structure beyond a simple definition of a void as 
an underdense region in the large-scale matter distribution. 
The tidal and the shear tensors encode more information 
than the density/velocity fields as they trace the collapsing 
or expanding nature of the matter field, which defines the 
dynamics of the Cosmic Web. 

The tracer that we use to define the voids is the frac¬ 
tional anisotropy (FA) computed from the set of eigenvalues 
of the tensor under consideration. The FA was initially in¬ 


troduced by Basser (1995 I to quantify the anisotropy degree 


of the diffusivity of water molecules through cerebral tissue 
in nuclear magnetic resonance imaging, thereby allowing to 
detect structures that restrict the otherwise isotropic Brow¬ 
nian movement of water molecules. Taking into account that 
this process is identical to the velocity shear in the large- 


scale matter distribution, Libeskind et al. (2013) introduced 


the concept of the FA in the context of Cosmic Web classi¬ 
fication schemes. 

In the next sections (|2]and © we establish the FA as 
a void tracer and then proceed to identify individual voids 
as basins of FA local minima. At this point we implement a 
watershed transform algorithm (Boucher & Lantuejoul 1979 


Beucher fc Meyer|1993 1 which has been used to define voids 


as catching basins of local minima of the density field ( |Platen| 
et al.|2007 Neyrinck|2008 1 . Finally we find and charachter- 
ize voids in a N-body m and © using density, velocity 
and fractional anisotropy profiles. We use these results to 
comment on the qualities of our algorithm (Sj 6 ]). 


2 ALGORITHMS TO FIND THE COSMIC 
WEB 

Our new void finding method is based on two existing cos¬ 
mic web classification schemes that work on cosmological 
N-body simulations. Both schemes depend on the construc¬ 
tion of tensors based on the Hessian of the potential (T- 
Web scheme) and the shear of the velocity (V-Web scheme). 
These algorithms have been used to develop other kind of 
studies such as the alignment of the shape, spin and pe¬ 
culiar velocity of dark matter halos with the cosmic web 
( Libeskind et al.|2013 Forero-Romero et al.|20i4 i. Here we 
summarize the most relevant aspects of each scheme. We re¬ 
fer the reader to the papers of Forero-Romero et al.| ( |2009| ) 
(T-Web) and Hoffman et al. 1 2012| ) (V-Web) for detailed 
descriptions. 


2.1 The tidal web (T-Web) 


This scheme was initially proposed by Hahn et al. (2007) 


to find the Cosmic Web based on the tidal tensor. The tidal 
tensor allows a classification in terms of the orbital dynamics 
of the matter field. This approach extends to second-order 
the equations of motion around local minima of the gravita¬ 
tional potential. The second-order term corresponds to the 
tidal tensor, which is defined as the Hessian matrix of the 
normalized gravitational potential 


T a p = 


d 2 (j) 

dx a dxp ’ 


(1) 


where the physical gravitational potential has been rescaled 
by a factor of AivGp in such a way that (f> satisfies the fol¬ 
lowing Poisson equation 


= S, 


( 2 ) 


with p the average density in the Universe, G the gravita¬ 
tional constant and 5 the dimensionless matter overdensity. 

Since the tidal tensor can be represented by a real and 
symmetric 3x3 matrix, it is always possible to diagonal¬ 
ize it and obtain three real eigenvalues Ai ^ A 2 ^ A 3 with 
its corresponding eigenvectors ui, uo, 113. The eigenvalues 
are indicators of the local orbital stability in each direction 
u;. The sign of the eigenvalues can be used to classify the 
Cosmic Web. The number of positive (stable) or negative 
(unstable) eigenvalues allows to label a location into one of 
the next four types of environment: voids (3 negative eigen¬ 
values), sheets ( 2 ), filaments ( 1 ) and knots ( 0 ). 

A modification to this scheme was introduced by |Forero^| 
Romero et al. (2009) by means of a relaxation of the stabil¬ 


ity criterion. The relative strength of each eigenvalue is no 
longer defined by the sign, but instead by a threshold value 
Ath that can be tuned in such a way that the visual impres¬ 
sion of the web-like matter distribution is reproduced. 


2.2 The velocity web (V-Web) 

The V-web scheme for environment finding introduced by 
Hoffma n et al.| |2012) is based on the local velocity shear 
tensor calculated from the smoothed dark matter velocity 
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Figure 1. Different fields in a slice of lh~ 1 Mpc width in the Bolshoi simulation. Top/bottom panels correspond to the T-Web/ V- Web. 
Fractional Anisotropy field (left); cosmic web classification (middle) for a given value of the threshold eigenvalue A t h where voids are 
white, sheets are light gray, filaments are dark gray and knots are black; individual voids (right) found by the a watershed algorithm on 
the Fractional Anisotropy field. Colours in the right panel correspond with the effective radius of each void. 


field in the simulation. This tensor is given by the following 
expression 


1 / dv a dvp A 

2H 0 \dxp dx a ) 


(3) 


where v a and x a represent the a component of the comoving 
velocity and position, respectively. Like the tidal tensor, T, a p 
can be represented by a 3 x 3 symmetric matrix with real 
values, making it possible to find three real eigenvalues and 
its corresponding eigenvectors. 

In this case we also use the relative strength of the three 
eigenvalues with respect to a threshold value A th to classify 
the cosmic web in the four web types already mentioned. 

Usually, the threshold A th is a free parameter that is 
tuned to reproduce the visual appearance of the comic web. 
However, in this paper this threshold does not play any role 
in our computations. In Figure [2] we offer a threshold es¬ 
timate (from the FA=0.95 as a function of Ai, this choice 
is explained in the next section) only as a guide to readers 
familiar with its meaning. 


3 A NEW VOID FINDING TECHNIQUE 


3.1 The fractional anisotropy 


The fractional anisotropy (FA), as developed by Basser 
(1995), was conceived to quantify the anisotropy of the dif- 


fusivity of water molecules through cerebral tissue in nuclear 
magnetic resonance imaging. In this manner, barriers such 
as microtubules and cell membranes, that restrict the other¬ 
wise isotropic Brownian motion of water molecules, can be 
detected. In a cosmological context, the problem of finding 
structures is very similar, where filaments and sheets re¬ 
strict the otherwise isotropic collapsing/expanding dynam¬ 
ics of the matter. Taking this into account, we present here 
the FA, much in the same way as Libeskind et al. ( |2013 1, to 
use it as a tracer of cosmic voids. 

The FA is defined as follows 


p » 1 /(Al - A 3 ) 2 + (A 2 - A 3 ) 2 + (Al - A 2 ) 2 

FA vJV-AfTIITAl-• (4 > 

where the eigenvalues can be taken from either the T-Web or 
the V-Web (FA-T-Web and FA-V-Web respectively). Such 
as it is defined, FA= 0 corresponds to an isotropic distri- 
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Figure 2. Fractional Anisotropy as a function of the overdensity (left) and the eigenvalue Ai (middle), and density contrast (right) as a 
function of the eigenvalue Ai for each web scheme (T-Web, continuous lines. V-Web, dashed lines). Fractional Anisotropy as a function of 
Ai (left), overdensity as a function of Ai (middle) and Fractional Anisotropy as a function of the overdensity (right). Thick central lines 
are the median and filled regions include 50% of the cells in each bin. The left panel shows the tight correlation between FA and Aj in 
voids. In the middle and right panel we show the equivalent eigenvalue threshold for the T-Web and V-web corresponding to FA=0.95. 
The right panel confirms that underdense regions can be identified with Fractional Anisotropy values less than 0.95. For this reason we 
define voids to be composed by cells with Fractional Anisotropy equal or less than this value. The threshold A, |, on the eigenvalues does 
not influence our results. We show it as a guide to readers familiar with the order of magnitude expectation for this quantity. 




FA 

Figure 3. Volume of each environment per FA (dVi/dFA, i = 
v,s,f,k) normalized by the total volume per FA ( dV/dFA ). 
Top/bottom corresponds to the T-Web/V-Web. For a given value 
of the FA, each line represents the volume contribution of the re¬ 
spective environment to the total volume with the same FA value. 
Above FA=0.95, only sheets and filaments contribute to all the 
volume of the simulation within that range of FA. 


bution (Ai = A 2 = A 3 ) and FA= 1 to a highly anisotropic 
distribution. 

I 11 the left and middle panels of Fig. [I] we show the FA 
field and web classification for both web schemes over a slice 
of an N-body simulation (described in Section|4|. Comparing 
these two panels we see that voids and knots (white and 
black in the middle panel of Fig. 1) display low FA values at 
their centres, becoming gradually more anisotropic at outer 
regions. On the other hand the filamentary structure (grey 
in the middle panel of Fig. 1) is traced by middle to high 
FA values. This characteristic is the key to use the FA as a 
tracer of cosmic voids. 


3.2 Fractional anisotropy as a void tracer 

Voids are regions where A 3 ^ A 2 ^ Ai ^ A t h- This 
implies that a void is completely fixed by the relative 
strength of the Ai eigenvalue with respect to the thresh¬ 
old. As we increase/decrease the threshold value Ath, voids 
increase/decrease progressively through contours of increas¬ 
ing/decreasing Ai. Voids are thus characterized by low val¬ 
ues of both FA and Ai. 

In Fig. [2] we show that the these two values are corre¬ 
lated. The left panel shows the correlation between the FA 
and Ai and the middle panel between the overdensity S and 
Ai. The right panel shows the correlation between 5 and the 
FA. This shows that the density has a larger scatter at fixed 
Ai than the FA, yet large and small values of the density 
remain still associated to low values of the FA. 

At this point, the eigenvalue Ai and the FA appear to be 
equally potential candidates for tracing cosmic voids, how¬ 
ever, the eigenvalue Ai presents, along with the overdensity, 
an undesirable feature for a void tracer. As we move out 
from inner parts of voids to outer parts, Ai and <5 increase 
monotonically making unclear where to set the void bound¬ 
ary, while for low density regions (5 < 1) the FA traces the 
boundary where anisotropic collapsing regions begin to dom¬ 
inate, i.e. at the ridge reached around <5 = 0 as shown in the 
right panel of Fig. [2] According to Fig. [3] these anisotropic 
areas are mainly associated to sheets and, to a lesser extent, 
filaments. 

From Figs. [T] and [2] we conclude that the FA is a good 
tracer of voids as it is almost perfectly correlated with low 
values of Ai. We propose that voids should be composed 
completely by regions of FA< 0.95. If we increase the values 
of Ai from its minimum until it we reach FA= 0.95 in [2] we 
find that this correspond to critical values of = 0.265 
and Ai" = 0.175 for the T-Web and V-Web, respectively. 
This means that setting A th to either Af /A} automatically 
produces voids with all the cells FA< 0.95. The middle pan¬ 
els in Fig. [I] show the web classification for this choice of 
Xth, demonstrating that this FA level is a reasonable choice 
to define voids. 
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3.3 Defining voids with a watershed algorithm 

The previous Section shows that FA is a good void tracer, 
but how should we actually define the boundary of individ¬ 
ual voids? For this purpose, we use the watershed transform 
algorithm (Beucher & Lantuejoul 1979| Beucher & Meyer 


19931 to identify a void as the basin of FA local minima 


with boundaries of FA= 0.95. The advantage of this defini¬ 
tion is that it does not require any assumption on the shape 
and/or morphology of the tentative voids. 

However, there are two main differences in our approach 
with respect to other watershed implementations. First, the 
watershed technique commonly uses the density field instead 
of the FA field as we do here (Platen et al.|[2007 Neyrinck 


|2008[ ). Second, we estimate all relevant quantities on a Carte¬ 
sian mesh of fixed cell size, while other works use an adaptive 


Delaunay tessellation (Schaap & van de Weygaert 2000). 
Nevertheless, in our case the average number of particles 
per cell is large enough (512 particles/cell) ensuring that 
the lowest density regions are sampled with at least ~ 10 
particles (<5 ~ —0.98), which seems to avoid spurious results 
from a noisy density field sampling. 

The watershed algorithm also needs a threshold value 
to reduce spurious features and prevent void hierarchization. 
If the density field is used, a typical threshold is 5 = —0.8 
(Platen et al. 2007), which means that, given two neigh¬ 


bouring voids, if the mean value of the density across their 
common boundary is below that threshold, the two voids 
are merged. In our case we have to find a corresponding FA 
value to define this threshold. 

In the right panel of Fig. [2] we see how the FA cor¬ 
relates with the matter overdensity. A value of 5 = —0.8 
corresponds to a value of FA~ 0.65 regardless of the web¬ 
finding scheme. This is the value we choose to remove ridges 
between two watershed voids. The right column in Fig. [T] 
shows all the individual voids that have been identified us¬ 
ing the watershed algorithm on the FA field. 


4 NUMERICAL SIMULATION 

We use the Bolshoi simulation to test our void finding 
method. This simulation follows the non-linear evolution 
of a dark matter density field on a cubic volume of size 
250ft~~ 1 Mpc sampled with 2048 3 particles. The cosmologi¬ 
cal parameters in the simulation are H m = 0.27, Ha = 0.73, 
h = 0.70, n = 0.95 and as = 0.82 for the matter density, cos¬ 
mological constant, dimensionless Hubble parameter, spec¬ 
tral index of primordial density perturbations and normal¬ 
ization for the power spectrum, respectively. These values 
are consistent with the ninth year of data of the Wilkin¬ 
son Microwave Anisotropy Probe (WMAP) ( |Hinshaw et al. 
2013). For more detailed technical information about the 


simulation, see |Klypin et al.| ( |2011| ) . 

We use data for the cosmic web identification that is 
publicly available through the MultiDark database http: 
//www.multidark.org/MultiDark/ which is described in 


obtain the data. For details see 

Forero-Romero et al. 

(2009) 

(T-Web); Hoffman et al. (2012); 

Forero-Romero et al. 

(2014 

(V-Web). This data is based on a cloud-in-cell (CIC) 

inter- 


onto a grid of 256 3 cells, corresponding to a spatial resolu¬ 
tion of 0.98/i _1 Mpc per cell side. These fields are smoothed 
with a gaussian filter with a width of a = 0.98/! _1 Mpc. The 
tidal and shear tensors and corresponding eigenvalues are 
computed through finite-differences over the potential and 
velocity fields. 


5 RESULTS 

We limit our results to voids with effective radius larger than 
the smoothing length of the density field, i.e. ~ lh -1 Mpc. 
For smaller interpolation/smoothing scales the shot noise in 
the low density regions starts to be noticeable. With that 
choice we find a void volume filling fraction 54.88% and 
47.06% for the FA-Vweb and the FA-Tweb, respectively. 
For larger smoothing scales (> 5/i _1 Mpc) the density and 
shear fields tend to be Gaussian and therefore follow closely 
the analytical environment abundance predictions for Gaus¬ 
sian random fields ( Forero-Romero et al.|2009 Alonso et al. 
20151. In this paper we limit ourselves to a smoothing scale 
~ 1 h -1 Mpc. 

In the following subsections we describe the results for 
the size distribution and different radial profiles for all our 
samples. 


5.1 The void size distribution 

Void shapes exhibit a wide range of geometries. To define 
their size we use its equivalent spherical radius or effective 
radius, defined as = [3/(4-7r) V] 1 / 3 , with V the total vol¬ 
ume of the void computed from the individual grid cells 
assigned to the void. In Fig. [T] we show the void size distri¬ 
butions for the T-Web and the V-Web. 

We see that the void distribution for the T-Web is 
broadly consistent with the expectations from a two-barrier 


polation of the density and velocity fields of the simulation 


problem (Sheth & van de Weygaert 20041. The formation 
of large voids is limited by the void-in-void mechanism 
(first barrier) where large voids are constituted hierarchi¬ 
cally of smaller ones. In turn, the formation of small voids 
is damped by the void-in-cloud mechanism (second barrier), 
where nearby collapsing structures limit the abundance of 
small embedded voids. 

We also find that the V-Web scheme produces an over¬ 
abundance of small voids compared to T-Web. A large num¬ 
ber of these small voids are embedded in overdense regions. 
They are already visible in the middle panel of Fig.[l]as tiny 
white bubbles located inside sheets. The existence of these 
small voids can be explained by dynamics of shell crossing in 
collapsing sheets as discussed in Hoff man et al.| ( [20 12]). The 
main argument is that as matter collides into a pancake-like 
fashion crossing sheets encounter each other at the symme¬ 
try plane. This effectively gives a positive divergence in the 
velocity field, resulting in a void identification by the V-Web 
algorithm. 

Comparing the abundance of large voids in the two web 
schemes, we find that the largest voids in the V-Web have 
r Bff ~ 10/i _1 Mpc, while the T-Web scheme includes voids 
as large as r eS ~ 15h -1 Mpc. This suggests that large voids 
in the T-Web scheme have a velocity structure that splits 
them in the V-Web. This increased granularity in the veloc- 
























































6 S. Bustamante and J.E. Forero-Romero 


Comoving volume [log 10 ( Mpc h 3 ] 



Effective radius fMocftl 


Comoving volume [log 10 ( Mpc h 1 ) ,! ] 

0.0 2.1 3.0 3.6 3.9 4.2 



Figure 4. Left panel, void size distribution. Right panel, fraction of overcompensated/subcompensated (C > 1/C < l)voids as a function 
of the effective radius. The continuous (dashed) curves correspond to the T-Web (V-Web). The abundance for small (r B ff < 2.0/r 1 Mpc.) 
overcompensated voids shoots up in the V-Web scheme. We provide in the body of the paper two reasons for this behaviour. 


ity structure of voids is evident in the left and right panels 
of Fig. □ 

5.2 Subcompensated and overcompensated voids 

We find that voids are distributed in two different families 
differentiated by the presence/absence of an overdense mat¬ 
ter ridge in their density profiles. To discriminate each void 
in one of the families we use the compensation index C It is 
defined as the mass of a void enclosed in a spherical volume 
of radius R and normalized by the mass of the same volume 
assuming it is filled by matter with the mean background 
density. 

c= w = 2 hi [,5(r)+1]r2dr (5) 

We choose an integration radius of R = 4r Bff , that is 
large enough to enclose the compensation ridge for a typical 
void in case there is one. This leads us to voids with C > 1 
having more mass than expected, constituting the family 
of overcompensated voids. These voids generally exhibit a 
compensation ridge associated to dense nearby structures. 
In the same fashion, voids with C < 1 constitute the family 
of subcompensated voids. 

In Fig. [5] we show the density and velocity profiles of 
voids split in these two families. In the left column it be¬ 
comes clear the difference between sub- and overcompen¬ 
sated voids. Tabic |7] lists the number of voids in each effec¬ 
tive radius bin for the two web finding methods. 

5.3 Density profiles 

We calculate the contrast density, radial-projected velocity 
and FA profiles. For this purpose we stack voids with similar 
effective radius. We also compute the profiles separately for 
subcompensated and overcompensated voids. For each void, 


Tweb Vweb 


r e ff [h 1 Mpc] 

C < 1 

C > 1 

C < 1 

C > 1 

< 2.0 

997 

495 

13480 

711 

2.0-3.2 

791 

471 

2581 

583 

3.2 - 4.4 

941 

469 

2168 

656 

4.4- 5.6 

869 

467 

1934 

724 

5.6 -6.8 

769 

395 

1287 

584 

6.8-8.3 

665 

362 

700 

315 

8.3- 12.0 

587 

333 

123 

61 

> 12.0 

29 

23 

0 

0 

total 

5648 

3015 

22273 

3634 


Table 1: Number of voids for each bin of effective radius and 
for subcompensanted and overcompansated samples. 


we take the distance of each member cell to the void centre 
along with the properties of interest. We normalize the radial 
distance with the effective radius and add the result to the 
stack. 

Fig- U shows the results of stacked density profiles for 
different void sizes. We calculate the profile out to a radius 
8 r.tf to capture the point where the overdensity reaches the 
mean value. 

The first interesting result is the overdensity value at 
the void’s centre. We find that larger voids have a lower over¬ 
density value. The largest voids (8.3 — 12 /i _1 Mpc) have an 
underdensity 5 ~ —0.95 while smaller voids (2—3.2 h -1 Mpc) 
fall around <5 ~ —0.8 at their centres. This holds for both web 
schemes. These values are consistent with most of the void 
finding schemes based on smooth and continuous fields from 


simulation or reconstruction procedures on surveys jPl ionis| 


& Basilakos|2002; Colberg et al. 2005 

Shandarin et al. 2006; 

Platen et al. 

2007| 

Neyrinck 2008| 

Munoz-Cuartas et al. 

2011 

ICeccare 

li et a 

. 2013 Paz et a 

. 2013; Neyrinck et al. 

2013 

Ricciardelli et al.120131, unlike geometrical approaches 
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Figure 5. Spherically averaged radial density profiles for stacked voids with different effective radii. The sample is split into subcom¬ 
pensated and overcompensated voids (right and left) for the web schemes T-Web and V-Web (top and bottom). 


based on point distributions, where central density values 
are generally higher (Colberg et al. [2008). 


A second feature about these profiles is their steepness 
at inner regions. In subcompensated voids, larger voids are 
steeper. Smaller voids exhibit moderate slopes, reaching the 
mean density at larger radii than larger voids. This sug¬ 
gests that small subcompensated voids are embedded into 
low density structures, while large subcompensated voids are 
surrounded by dense structures, reaching the mean density 
at lower effective radii than smaller voids. In the overcom¬ 
pensated case larger voids reach first both the compensation 
ridge and then the mean density value. 


Regarding overcompensated voids, a final result is re¬ 
lated to the height of the compensation ridge: the larger 
the void size, the lower the ridge height. This implies that 
overcompensated smaller voids are embedded in very high 
density regions, unlike their subcompensated counterpart, 
thus indicating two possibly different processes for small 


voids formation. Larger voids exhibit lower ridges as outer 
radial layers also includes all sort of structures, thus being 
the difference between large overcompensated and subcom¬ 
pensated voids less conclusive. 


All the previous results hold for both finding schemes. 
This suggests an universal behaviour for the radial den¬ 
sity profile in two families of subcompensated and over- 
compensanted voids. This goes in the same direction of re- 


rdelli et al. 

2013 

1 and external structure of voids 1 

Lavaux 

& Wandelt 

2012 

|Ceccarelli et al. 2013; Paz et a 

1. 20131 


Hamaus et al. 20141. Indeed, our results extend the find- 
ings of |Hamaus~et al. (2014j) into the range of voids with 
size r aff < 10h“ i Mpc. 
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Normalized radius r/r eff 


Normalized radius r/r cff 


Figure 6. Spherically averaged radial velocity profiles for stacked voids with different effective radii. The sample is split into subcom¬ 
pensated and overcompensated voids (right and left) for the web schemes T-Web and V-Web (top and bottom). 


5.4 Velocity profiles 

In Fig. [6] we present the radial velocity profiles. Positive 
values correspond to outflows with respect to the centre. 

We find that subcompensated voids have outflowing ve¬ 
locity profiles all the way up to the radius where the aver¬ 
age radial density reaches zero. In voids with sizes r e s < 
8/i -1 Mpc the outflow is always positive, consistent with the 
fact that their density profiles do not reach the level 5 = 0 
in the range of explored radii. This behaviour indicates that 
matter is being pulled out of the void into external higher 
density features. 

On the other hand, overcompensated voids initially ex¬ 
hibit outward profiles, (as expected from a low density re¬ 
gion) and approximately at the radius of the compensation 
ridge, the velocity reaches a peak, decreases and becomes 
negative, showing the infalling flow of matter further than 
the compensation ridge. This shows that the high density 


structures associated to the compensation ridge dominate 
the matter flow both from inside and outside the void. 

As in the density, these results are also consistent with a 
Universal velocity profile for voids ( |Paz et ah 2013i [H amaus| 
et al.|2014 i. 

5.5 Fractional Anisotropy profiles 

Fig-0 shows the results for the FA profiles. We find that 
the FA clearly magnifies the difference between the inter¬ 
nal profile r/r e ff < 1 and the external profile r/r eff > 1. 
Subcompensanted voids in the T-Web reach an asymptotic 
background value of F=0.8 almost right at r/r eff = 2. For 
the V-Web the results are similar, albeit with a slower trend 
with the effective radius. Overcompensated voids reach a 
maximum FA at the same effective radius r/r e tr = 1 and 
decline to reach an asymptotic value of FA=0.85 at larger 
radii. 
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Figure 7. Spherically averaged Fractional Anisotropy profiles for stacked voids with different effective radii. The sample is split into 
subcompensated and overcompensated voids (right and left) for the web schemes T-Web and V-Web (top and bottom). 


The central FA values also show a magnified trend with 
the void size. Large voids have lower FA values, amplify¬ 
ing the same trend observed with the central density value. 
Finally, we observe that voids in the V-Web scheme span a 
larger range of central FA values than in the T-Web scheme. 

The difference between the radius where the density 
ridge is reached (r/r e ff = 3) and the radius of the FA max¬ 
ima (r/r B ft = 1, most manifest in the T-Web voids) justi¬ 
fies in a more quantitative way the qualitative argument we 
present in sub-section |3.3| to define void boundaries in our 
method. Namely that as we increase the void boundary mov¬ 
ing away from the centre, we reach first middle density walls, 
before reaching high density structures associated with high 
anisotropy. This implies that voids defined with the FA never 
present positive overdensities, making the FA maximum a 
reasonable boundary for voids compared to the traditional 
definition that puts the boundary at the density ridge. 

We recognize that our choice produces smaller voids as 


compared with other voids finding methods. We consider 
that this has the advantage of avoiding the contamination 
from external structures. 


6 CONCLUSIONS 

In this paper we propose that the anisotropy of the eigenval¬ 
ues from tidal and velocity shear tensors is a good tracer of 
cosmic voids. Based on this idea we go on to implement a wa¬ 
tershed algorithm on the fractional anisotropy to find voids 
in a N-body cosmological simulation. We perform different 
test on the voids found on the anisotropy of two tensorial 
schemes, the T-Web and the V-web. 

The first quantity we examine is the void size distribu¬ 
tion characterized by an effective radius. We find that for the 
T-Web scheme the void size distribution has a shape consis¬ 
tent with standard expectations from a two-barrier setup. 
However, for the V-Web we find an overabundance of the 
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smallest voids. This is due to two reasons. First, there are 
artefacts in the web finding scheme inside sheets. Second, 
there V-web produces smaller fragmented voids inside the 
largest voids found in the T-Web classification. This void 
fragmentation hints at a complex velocity structure in large 
voids compared to its simpler density/gravitational poten¬ 
tial anatomy. 

A second step in the void characterization is the separa¬ 
tion into subcompensated and overcompensated samples. In 
the T-Web 35% of the voids are overcompensated and 65% 
subcompensated, meaning that they are located in denser 
regions with a clear delimiting ridge. For the V-web the sub¬ 
compensated fraction rises to 85% as the void size decreases, 
supporting the picture of smaller spurious voids located in 
medium density matter sheets. 

Finally we proceed with a more detailed void character¬ 
ization through the radially averaged profiles of density, ra¬ 
dial velocity and fractional anisotropy. In this case we study 
separately the subcompensated and overcompensated voids, 
split this time in samples with different effective radii. 

For the density the most interesting feature is that the 
profile has a similar shape for all the different void sizes 
once the radial coordinate is expressed in units of the ef¬ 
fective radius. In these profiles it is evident the presence of 
an overdense ridge around 2 < r/r B ft < 3 for the overcom¬ 
pensated voids. Subcompensated voids do not show such a 
ridge keeping its density lower than the average value up to 
r/r. t{ ~ 8. 

In other studies where voids are found using considera¬ 
tions on the density fields the ridge is, almost by definition, 
located around r/r e ft ~ 1 and the subcompensated voids 
reach average density around r/r e a ~ 3. This indicates that 
voids defined by FA boundaries are smaller than voids found 
by density only considerations close by a radial factor of 2. 

The velocity profiles show the expected correlation with 
the features already observed in the density. Most notably 
the radial velocity goes to zero close to flat regions in the 
density profile and is positive around regions of increasing 
density with radius. 

The FA profiles also show remarkable similarities inside 
each class of subcompensated and overcompensated voids 
and a close correlation with the behaviour in the density 
and velocity profiles. From these profiles it becomes clear 
that the limit of overcompensated void coincides with a ridge 
in of FA values close to 0.9. This value is close to the limit 
FA=0.95 we impose in the watershed algorithm, indicating 
that these voids are very close to spherical. For the sub¬ 
compensated voids the FA quickly increases to FA values of 
~ 0.7 — 0.8 at r/r aff ~ 1 to hnd a plateau, which is quite far 
from the limiting FA= 0.95 suggesting that these voids are 
ellipsoidal. 

Put together, the results for the radially averaged pro¬ 
files support the evidence for a universal density profile. Our 
results extend to smaller void sizes t he work done by |Hamaus| 
|et al.| p014| on simulations and by |Ceccarelli et ah (20131 
on SDSS data. 

The method we present here finds voids with reasonable 
properties compared to different results published in the lit¬ 
erature. This gives support to a new tracer, the Fractional 
Anisotropy, to be used in the study of cosmic voids. Fur¬ 
thermore, the complementary physical picture in the two 
different web finding algorithms opens the possibilities to 


make a joint analysis of the tidal and dynamical structure 
of voids. 
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